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Abstract 


To  investigate  the  mechanical  and  other  condensed  phase  properties  of  energetic 
materials  using  atomistic  simulation  techniques,  the  COMPASS  force  field  has  been 
expanded  to  include  high-energy  nitro  functional  groups.  This  report  presents  the 
parameterization  and  validation  of  COMPASS  for  nitrate  esters  (-ONO2).  The  functional 
forms  of  this  force  field  are  of  the  consistent  force  field  type.  The  parameters  were 
derived  with  an  emphasis  on  the  nonbonded  parameters,  which  include  a  Lennard-Jones 
9-6  function  for  the  van  der  Waals  (vdW)  term  and  a  Coulombic  term  for  an  electrostatic 
interaction.  To  validate  the  force  field,  molecular  mechanics  calculations  and  molecular 
dynamics  simulations  have  been  made  on  a  variety  of  molecules  containing  the  nitrate 
ester  functionality.  Using  this  force  field,  excellent  agreement  has  been  obtained  between 
the  calculated  and  experimental  values  for  molecular  structures,  vibrational  firequencies, 
liquid  densities,  heats  of  vaporization,  crystal  stmcture,  mechanical  properties  and  lattice 
energy. 
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1.  Introduction 


The  ballistic  performance  of  a  propellant  in  a  gun  system  is  known  to  be  dependent  upon  the 
mechanical  response  of  the  propellant  grains  to  the  high  stress  and  high  strain  environment 
experienced  during  the  interior  ballistic  cycle.  In  addition,  although  not  totally  understood,  the 
vulnerability  of  a  propellant,  which  is  the  response  of  the  propellant  to  a  variety  of  stimulii  such 
as  impact  and  heat,  is  known  to  be  influenced  by  the  mechanical  properties  of  the  material. 
Generally  speaking,  poor  mechanical  properties  yield  poor  vulnerability  results;  however,  good 
mechanical  properties  don’t  always  yield  good  vulnerability  results.  In  numerical  models  of  the 
interior  ballistic  process,  the  propellant  grain  is  considered  to  be  an  incompressible, 
nondeformable  solid  that  bums  with  a  mass  generation  rate  that  is  determined  by  the  measured 
bum  rate  and  the  exposed  surface  area.  If  grain  fracture  occurs,  more  surface  area  is  exposed, 
which  results  in  an  increase  in  the  expected  mass  generation  rate.  This  can  lead  to  decreased 
performance,  or  in  the  worst  case  scenario,  complete  failure.  Our  objective  in  initiating  this 
study  is  to  develop  the  capability  to  use  molecular  modeling  techniques  to  characterize  the 
physical  properties  of  high-energy  materials  at  the  atomistic  level.  The  information  gained  from 
this  modeling  effort  can  help  direct  changes  in  the  formulation  and  the  processing  techniques 
used  in  the  manufacturing  of  the  propellant,  which  may  ultimately  result  in  improving  the 
performance  and  vulnerability  properties  of  the  propellant. 

For  modeling  the  properties  of  large  molecular  systems  and  condensed  phases,  molecular 
modeling  techniques  that  are  force  field  based  have  distinct  advantages  over  ab  initio 
methodologies.  This  is  not  only  because  the  force  field  method  is  several  orders  of  magnitude 
faster  than  any  ab  initio  method,  but  also  because  the  information  that  an  ab  initio  method 
provides  is  often  not  necessary  for  these  applications.  The  properties  of  interest  are  generally  the 
result  of  the  statistical  average  of  atomistic  movement  over  a  much  longer  time  scale  than  the 
rapid  electron  motion  that  an  ab  initio  method  describes.  In  addition,  the  most  important 
interaction  terms  in  simulating  the  condensed  phase  are  the  nonbonding  forces  (in  particular,  the 
dispersion  forces),  which  are  extremely  difficult  to  accurately  describe  using  ab  initio  methods. 
Reports  have  recently  appeared  which  document  the  accurate  prediction  of  a  wide  variety  of 
properties  of  molecules  in  both  the  condensed  phase  and  in  isolation  using  the  COMPASS 
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(condensed-phase  optimized  molecular  potentials  for  atomistic  simulation  studies)  force  field 
[1—2].  Popular  force  fields,  such  as  MM3  [3-13],  AMBER  [14-15],  and  CHARMM  [16],  have 
been  designed  mainly  to  study  biologically  interesting  molecules.  The  COMPASS  force  field, 
on  the  other  hand,  has  been  specifically  designed  for  material  science  applications.  It  is  a  class  n 
ab  initio  force  field  in  that  it  employs  complex  fimctional  forms  and  is  derived  firom  extensive  ab 
initio  data.  Consequently,  it  can  be  used  to  accurately  predict  several  molecular  properties, 
including  molecular  structures,  conformations,  and  vibrations.  The  nonbonded  parameters  in 
COMPASS  have  been  optimized  using  condensed-phase  (liquid  and  crystal)  data  so  that  several 
thermophysical  properties  of  molecular  liquids  and  crystals  can  be  well  reproduced.  It  is 
important  to  note  that  COMPASS  has  been  developed  to  model  the  static  properties  of 
condensed  phase  molecules.  It  is  not  designed  to  study  reactions;  the  development  of  a  reactive 
force  field  is  a  significant  undertaking  and  is  outside  the  scope  of  this  work.  COMPASS  has 
broad  coverage  in  the  major  application  areas  of  material  science.  It  has  been  parameterized  to 
study  most  common  organic  molecules,  organic  and  inorganic  polymers,  zeohtes,  and 
metal/transition-metal  oxides.  However,  some  of  the  fimctional  groups  required  to  model 
energetic  materials  have  not  until  now  been  parameterized  and  included  in  the  COMPASS  force 
field.  For  instance,  a  common  propellant  is  composed  of  three  principal  energetic  ingredients: 
nitrocellulose  (NC,  C6H7.55N2.45O9.90),  nitroglycerin  (NG,  O2NOCH2-CHONO2-CH2ONO2),  and 
diethyleneglycol  dinitrate  (DEGDN,  02N0-(CH2)2-0-(CH2)2-0N02).  In  the  propellant  blend, 
the  NC  is  present  as  an  energetic  binder  and  NG  serves  as  an  energetic  gelling  agent.  DEGDN  is 
also  an  energetic  ingredient;  however,  it  also  plays  the  role  of  a  plasticizing  agent  in  the 
propellant.  The  common  chemical  feature  in  each  of  these  compounds  is  the  nitrate 
ester  (-ONO2)  fimctional  group.  Compounds  containing  nitro  (-NO2)  are  currently  parameterized 
in  the  COMPASS  force  field  and  contain  most,  but  not  all,  of  the  parameters  needed  to  model  the 
nitrate  esters.  This  work  reports  on  the  parameterization  and  validation  of  the  additional  atom 
types  needed  to  model  the  nitrate  ester  fimctional  group  in  the  COMPASS  force  field. 

2.  Force  Field  Development 

2.1  Ab  Initio  Calculations.  Very  few  high  quality  measurements  of  the  molecular  structures 
and  conformations  have  been  reported  for  high  energy  compounds  containing  the  nitrate  ester 
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functional  group.  Since  experimental  data  is  lacking,  the  validation  of  the  force  field  has  to  be 
based  on  accurate  ab  initio  data.  In  order  to  determine  the  level  of  calculation  that  is  reasonable 
both  in  terms  of  accuracy  and  computational  expense,  we  carried  out  an  evaluation  of  several 
candidate  ab  initio  methods.  This  was  done  by  performing  a  series  of  calculations  on  two  model 
nitrate  ester  compounds,  methyl  and  ethyl  nitrate,  as  shown  in  Figure  1,  and  included  both  the 
trans  and  gauche  conformers  of  ethyl  nitrate.  The  following  levels  of  theory  were  used: 
HF/6-31g(d),  HF/6-311g(d,p),  MP2/6-311g(d,  p),  MP2/6-311g(2df,  2p),  SVWN/6-311g(d,  p), 
and  B3LYP/6-311g(d,  p).  These  calculations  were  carried  out  using  the  GAUSSIAN  94  (G94) 
[17]  suite  of  quantum  chemistry  programs  running  on  a  Silicon  Graphics  Origin  2000  located  at 
the  U.S.  Army  Research  Laboratory  Major  Shared  Resource  Center.  The  calculations  were 
completed  using  the  default  convergence  criteria  for  all  methods  and  the  default  grid  sizes  for  the 
DFT  methods  contained  in  G94.  Table  1  lists  the  calculated  results  at  the  different  levels  of 
theory  and  the  experimental  data  for  methyl  nitrate,  while  Tables  2  and  3  list  the  calculated 
results  and  experimental  data  for  the  trans  and  gauche  conformers  of  ethyl  nitrate,  respectively. 
The  calculated  energies  of  the  ethyl  nitrate  conformers  are  presented  in  Tables  2  and  3.  In 
addition,  the  energy  difference  between  the  trans  and  gauche  conformer  at  each  level  of  theory  is 
presented  in  Table  3. 

A  comparison  of  the  calculated  structure  of  methyl  nitrate  at  the  various  levels  of  theory  with 
the  experimental  stmcture  determined  by  Cox  and  Waring  [18]  yields  the  following  results:  the 
HF/6-31g(d)  and  HF/6-31  lg(d,p)  bond  lengths  are  within  2%  of  experiment,  and  the  bond  angles 
are  within  0.9%  of  experiment;  the  MP2/6-311g(d,p)  and  the  MP2/6-311g(2df,2p)  bond  lengths 
are  within  0.4%  of  experiment,  and  the  bond  angles  are  within  0.4%  of  experiment;  the 
SVWN/6-31  lg(d,p)  bond  lengths  are  within  1.2%  of  experiment,  and  the  bond  angles  are  within 
0.8%;  and  finally,  the  B3LYP/6-31  lg(d,p)  bond  lengths  and  bond  angles  are  both  within  0.5%  of 
experiment.  Similar  agreement  is  observed  between  our  calculations  and  the  results  obtained 
fi'om  the  microwave  studies  of  Scroggin  et  al.  [19]  on  ethyl  nitrate.  An  analysis  of  the  calculated 
energies  of  ethyl  nitrate  shows  the  difference  between  the  isomers  to  be  0.73  kcal/mole  at 
HF/6-31g(d),  0.76  kcal/mole  at  HF/6-31  lg(d,p),  0.18  kcaVmole  at  MP2/6-311g(d,p), 
0.22  kcal/mole  at  B3LYP/6-311g(d,p),  and  0.13  kcal/mole  at  MP2/6-311g(2df,2p).  In  each  of 
these  calculations,  the  trans  conformer  is  predicted  to  be  more  stable  than  the  gauche  conformer. 
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Figure  1.  Methyl  and  Ethyl  Nitrate. 
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Table  1.  Optimized  Geometries  of  Methyl  Nitrate 


2.  Optimized  Geometries  of  the  Trans  Conformer  of  Ethyl  Nitrate 


Table  3.  Optimized  Geometries  of  the  Gauche  Conformer  of  Ethyl  Nitrate 


These  results  are  consistent  with  the  calculations  reported  by  Durig  and  Sheehan  [20].  The 
exception  in  our  series  of  calculations  is  at  the  SVWN/6-311g(d,p)  level  of  theory  where  the 
gauche  isomer  is  calculated  to  be  the  more  stable  of  the  two  by  0.6  kcal/mole.  The 
B3LYP/6-31  lg(d,p)  method  yields  good  results  for  the  molecular  structures  of  the  model  nitrates 
and  for  the  conformational  energy  differences  of  ethyl  nitrate,  with  considerably  lower 
computational  requirements  than  the  MP2  calculations.  This  method  provides  an  ideal 
combination  of  accuracy  and  efficiency,  consequently,  the  B3LYP/6-311g(d,p)  method  was  used 
throu^out  the  remainder  of  this  study  for  validating  the  structures  of  additional  nitrate  esters, 
such  as  NC,  NG,  and  DEGDN,  for  which  little  experimental  data  is  available. 


2.2  Parameterization  Method.  The  fimctional  forms  of  the  COMPASS  1 .0  force  field  are 
the  same  as  the  CFF-type  force  fields  [21-29],  as  in 
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The  potential  flmctions  can  be  divided  into  two  categories — ^valence  terms  including  the 
diagonal  and  off-diagonal  cross  coupling  terms  and  the  nonbonded  interaction  terms.  The 
valence  terms  include  E^,  E^  and  E^  for  bond,  angle,  torsion,  and  out-of-plane  angle 
coordinates,  respectively,  and  Ef^y,  Ef^g,  Ef^^,  Egg,,  zadEgg^  represent  the  cross-coupling  terms 

between  internal  coordinates.  The  cross-coupling  terms  are  important  for  predicting  vibrational 
frequencies  and  structural  variations  associated  with  conformational  changes.  Among  the  cross¬ 
coupling  terms,  the  bond-bond  Ej^g,  bond-angle  Egg,  and  bond— torsion  Eg^  are  the  most 

significant.  The  nonbonded  terms,  which  include  a  “soft”  Leimard-Jones  9-6  (L-J)  potential  for 
the  van  der  Waals  (vdW)  interaction  and  a  Coulombic  term  for  the  electrostatic  interactions,  are 
used  for  interactions  between  pairs  of  atoms  that  are  separated  by  three  or  more  intervening 
atoms,  or  those  that  belong  to  different  molecules.  The  L-J  parameters  (e  and  r°)  for  like  atom 
pairs  are  adjustable  parameters.  For  unlike  atom  pairs,  a  sixth  order  combination  law  [30]  is 
used  to  calculate  the  off-diagonal  parameters  as  follows: 


and 


(2) 


(3) 


The  electrostatic  interaction  is  represented  by  the  partial  atomic  charge  model  using  charge 
bond-increments,  dy,  as  a  measxue  of  the  charge  separation  between  two  valence-bonded  atoms. 
The  net  partial  charge  of  an  atom,  is  obtained  as  a  summation  of  all  charge-bond  increments 
related  to  this  atom: 


9/ 


J 


(4) 
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The  details  on  the  methodology  of  parameterization  are  reported  elsewhere  [1,2].  Basically, 
the  valence  parameters  (both  diagonal  and  off-diagonal  cross-coupling  terms)  and  charge-bond 
increments  were  derived  by  least-squares  fitting  to  the  HF/6-31g(d)  data  calculated  for  the  model 
compounds  methyl  and  ethyl  nitrate.  The  ab  initio  data  includes  electrostatic  potentials, 
energies,  and  the  first  and  second  derivatives  of  the  energies.  The  resulting  parameters  were 
subsequently  scaled  by  a  set  of  generic  factors  to  correct  the  systematic  errors  of  the  HF/6-31g(d) 
calculations.  The  resulting  force  field  (quantmn  mechanics  force  field  or  QMFF)  was  then 
systematically  validated  and  modified  to  fit  experimental  or  higher-level  ab  initio  data  for 
molecules  in  isolation.  The  vdW  nonbonded  parameters  (L-J  9-6  terms)  were  initially 
transferred  firom  other  organic  systems  [21,  22,  25-28,  31,  32].  After  the  valence  parameters 
were  derived,  they  were  subject  to  optimization  using  MD  simulations  of  liquids.  The  whole 
validation  and  optimization  procedure  was  repeated  until  a  consistent  fit  was  obtained  for  both 
the  gaseous  and  condensed  phases. 

2.3  Molecular  Dynamics  Simulations.  Molecular  dynamics  (MD)  simulations  were  carried 
out  using  the  software  package  hisightll/Discovery.  For  liquids  and  crystals,  a  periodic  cell  with 
explicit  miniTmim  image  convention  [33,  34]  was  built  for  each  of  the  model  compounds  studied. 
The  cubic  cell  edges  ranged  in  length  firom  20  to  30  A  and  contained  1000  to  1500  atoms.  A 
charge-group-based  cutoff  method  with  tail  correction  was  used  to  evaluate  the  nonbonded 
interactions  in  all  of  the  liquid  simulations.  It  is  assumed  in  the  charge-group-based  cutoff 
method  that  the  radial  distribution  functions  converged  to  unity  beyond  the  cutoff  distance 
[33,  34],  which  in  our  simulations  was  9.5  A.  For  the  crystal  simulations,  the  Ewald  summation 
method  [33,  34]  was  used  for  both  the  vdW  and  the  electrostatic  terms.  Constant  volume  and 
temperature  (NVT)  ensembles,  with  a  velocity  Verlet  integrator  [33, 34]  and  Andersen’s  [35, 36] 
temperature  control  method,  were  used  for  the  parameterization.  Constant  pressure  and 
temperature  (NPT)  simulations  were  carried  out  using  a  modified  velocity  Verlet  [33,  34] 
integrator  with  the  Berendsen  [37]  pressure  control  method  for  the  validation  calculations.  A 
time  step  of  1  fs  was  used  in  all  of  the  MD  simulations. 

The  initial  configuration  of  the  liquid  being  simulated  was  constructed  using  the  following 
procedure.  First,  molecules  were  uniformly  placed  into  a  cubic  cell  that  was  large  enough  so  that 
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there  was  no  strong  repulsion  between  any  two  molecules.  The  system  was  then  randomized  by 
running  a  NVT  simulation  for  several  hundred  steps  at  a  temperature  of  2000  K,  which  is  well 
over  the  boiling  point  of  all  of  the  liquids  studied.  Following  randomization,  the  system  was 
gradually  compressed  to  the  target  density  in  a  high  pressure  (5,000-50,000  bar)  NPT 
simulation.  Finally,  a  pre-equilibrium  process  was  performed  using  a  simulated  annealing 
technique,  during  which  the  temperature  was  gradually  reduced  from  2000  K  to  300  K.  The  pre¬ 
equilibration  took  about  50-100  ps,  which  is  usually  adequate  for  liquids  of  small  molecules 
[1,2].  The  average  periods  were  50  ps  for  NVT  simulations  and  100  ps  for  NPT  simulations. 

3.  Parameterization 

3.1  Valence  Parameters.  As  stated  previously,  the  valence  parameters  were  derived  from 
ab  initio  data  using  a  least  squares  fitting  procedure  [1,  2].  Simple  nitro-containing  (-NO2) 
compounds  (such  as  nitrobenzene  and  nitro-alkanes),  which  share  most  of  the  same  parameters 
with  the  nitrate  esters,  were  parameterized  in  the  COMPASS  1.0  force  field  [1,  2].  In  this 
project,  additional  model  compounds,  methyl  and  ethyl  nitrate,  were  calculated  at  the  HF/6- 
31g(d)  level  of  theory  to  sample  the  “missing”  interaction  energy  terms.  To  maintain 
consistency  with  the  other  COMPASS  parameters,  the  HF/6-31g(d)  data  were  used  to  derive  the 
missing  valence  parameters,  which  due  to  historic  limitation,  were  derived  at  the  same  level  of 
theory  since  the  early  1990’s.  The  systematic  errors  between  the  HF/6-31g(d)  calculations  and 
experimental  measurements  have  been  well  documented  [21-29],  and  a  set  of  scaling  factors 
that  correct  the  errors  have  been  established  and  implemented  throughout  the  development  of 
COMPASS.  With  the  additional  data,  the  remaining  valence  parameters  were  derived  with  an 
accuracy  level  similar  to  the  other  systems  that  have  been  parameterized  and  reported  previously 
[1,2].  The  new  parameters  were  then  subjected  to  the  same  empirical  adjustment  procedure  to 
yield  a  better  overall  fit  of  the  stmctural  and  conformational  properties  of  the  molecules  in 
isolation.  The  final  parameters  for  nitro-containing  compounds  are  given  in  the  Appendix. 
Three  atom  types  [1,  2]  are  used  for  these  compoimds:  n3o  for  the  nitrogen  atom  of  the  nitro 
group,  0I2  for  the  oxygen  atoms  in  the  nitro  group,  and  o2n  for  the  ester  oxygen  in  the  nitrate 
ester.  Note  that  only  n3o  is  a  specific  atom  type  for  the  nitro  moiety;  both  ol2  and  o2n  are 
equivalent  to  the  more  generic  atom  types,  o2  and  ol,  for  the  valence  terms  [1,2].  Together  with 
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the  atom  types  and  parameters  defined  for  alkanes  and  benzenes  [1,  2],  the  parameters  provide 
complete  coverage  of  simple  nitro-containing  compounds  and  nitrate  esters. 

3.2  Atomic  Partial  Charges.  As  a  significant  part  of  the  parameterization  of  the  nonbonded 
terms,  atomic  partial  charges  have  drawn  considerable  attention  during  the  development  of  the 
COMPASS  force  field.  Generally  speaking,  a  constrained  fitting  of  the  electrostatic  potential 
[1,2]  (CESP)  is  used  to  derive  the  atomic  partial  charges.  This  approach  ensures  the 
transferability  of  the  charge  parameters,  while  maintaining  an  overall  best  fit  of  the  ab  initio 
electrostatic  potential  energy  surfaces.  The  actual  partial  charges  for  nitrous  acid,  methyl  nitrite, 
nitrobenzene,  methyl  nitrate,  and  ethyl  nitrate,  obtained  irom  both  Mulliken  analysis  and  CESP 
fitting,  are  listed  in  Table  4  for  comparison.  The  CESP  charges  are  represented  by  bond  charge 
increments  (6,y)  which  represent  the  charge  flow  between  two  adjacent  atoms,  i  and  j,  as  given  in 
the  Appendix.  Several  parameters  (for  alkanes  and  benzenes)  are  transferred  and  fixed  during 
the  CESP  fitting  procedure;  a  single  bond-charge  increment  applies  to  all  compoimds  containing 
the  same  chemical  bond  (eg.,  one  b„3o,oi=  for  all  N-O  bonds  in  the  nitro  group  and  the  nitrate 
group).  Consequently,  the  partial  charges  obtained  are  more  “symmetric”  than  those  obtained 
from  the  Mulliken  analysis.  Despite  the  enforced  constraints  in  the  CESP  fit,  the  overall  quality 
of  the  reproduced  ab  initio  electrostatic  potential  energy  surfaces  is  satisfactory.  This  can  be 
measured  by  calculating  the  root  mean  square  (RMS)  errors  (in  kcal/mol)  between  the  ab  initio 
electrostatic  potentials  and  the  potentials  calculated  using  the  partial  charges.  The  RMS  errors 
are  presented  in  Table  5.  Using  the  CESP  charges,  the  calculated  RMS  errors  in  the  ESP  are 
significantly  smaller  than  those  calculated  using  the  Mulliken  charges.  An  additional  evaluation 
of  the  electrostatic  parameters  is  given  in  Table  6,  where  the  molecular  dipole  moments 
calculated  using  the  Mulliken  and  CESP  charges  are  compared  against  experimental  and  ab 
initio  results.  Although  the  dipole  moments  were  not  used  as  input  data  for  fitting,  the  CESP 
charges  reproduce  these  values  reasonably  well. 

3.3  The  van  der  Waals  (vdW)  Parameters.  With  the  valence  and  charge  parameters  fully 
determined,  the  only  remaining  terms  to  be  addressed  are  the  weak  vdW  interactions.  These 
terms,  which  are  known  to  play  a  critical  role  in  the  simulation  of  condensed  phases,  were 
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Table  4.  Comparison  of  Atomic  Partial  Charges 
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Table  5.  RMS  Errors  (kcal/mol)  Between  Ab  Initio  ESPs  and  ESPs  From  Partial  Charges 


Molecule 

Mulliken 

CESP 

HNO2 

2.397 

0.735 

CH3NO2 

3.491 

0.771 

C6H5NO2 

3.666 

1.367 

CH3ONO2 

2.502 

0.860 

C2H5ONO2 

2.840 

0.951 

Table  6.  Comparison  of  Calculated  and  Experimental  Molecular  Dipole  Moments 


Molecule 

HF/6-31g(d) 

Mulliken 

HNO2 

3.028 

3.978 

3.038 

— 

CH3NO2 

4.048 

5.958 

3.994 

— 

C6H5NO2 

5.031 

7.635 

4.574 

— 

CH3ONO2 

3.835 

5.060 

3.896 

3.081“ 

C2H5ONO2 

4.125 

5.374 

3.931 

3.39”  (trans) 
3.23’’  (gauche) 

“Reference  [18]. 
’deference  [19], 


subject  to  parameterization  using  condensed  phase  data.  The  Lennard- Jones  9-6  (L-J)  function 
is  used  for  representing  the  van  der  Waals  energies.  Most  of  the  noribonded  vdW  parameters  are 
transferred  from  alkanes,  ethers,  and  nitro  compounds  [1,2].  In  this  project,  only  one  new  atom 
type  (o2n)  was  introduced,  which  requires  two  new  vdW  parameters  (the  L-J  well  depth,  e,  and 
the  vdW  radii,  i")  to  be  defined.  We  selected  liquid  methyl,  ethyl,  propyl,  isopropyl,  and  butyl 
nitrate  to  be  used  in  the  determination  of  these  parameters.  The  thermophysical  properties  of 
density,  p,  and  the  heat  of  vaporization,  AHv,  of  these  liquids  were  calculated.  The  heat  of 
vaporization  was  obtained  from  the  cohesive  energy  density  (Eced)  by  use  of  the  followmg 
equation: 
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MJy  =Eced-MI  p  +  RT , 


(5) 


where  M  is  the  molecular  wei^t  and  p  is  the  density.  The  parameterization  process  consisted  of 
repeated  MD  simulations  during  which  the  adjustable  parameters  (e  and  r®)  were  modified.  By 
comparing  the  calculated  pressures  and  the  cohesive  energies  at  a  given  temperature  against  the 
experimental  data,  the  optimized  parameters  were  identified  numerically  fi'om  5  to  10  data 
points. 

4.  Validation 

4.1  Gas  Phase  Molecular  Properties.  The  validation  of  molecular  properties  for  molecules 
in  the  gas  phase  was  based  on  the  molecular  mechanics  calculations.  These  calculations  were 
performed  on  the  isolated  (i.e.,  gas  phase)  molecules  of  methyl  and  ethyl  nitrate.  In  each  case, 
the  calculations  consisted  of  full-energy  minimizations  followed  by  calculation  of  the  Hessian 
matrix.  These  results  were  compared  with  either  the  experimental  data  or  with  the  results  firom 
high-level  ab  initio  calculations. 

4A.1  Structures.  The  most  basic  property  to  predict  is  the  structure  of  the  molecule.  It  is 
well  known,  for  example,  that  a  small  deviation  in  the  bond  length  can  have  a  potentially 
significant  effect  on  the  liquid  or  crystal  density  obtained  fi:om  an  MD  simulation. 
Consequently,  prior  to  running  the  condensed  phase  simulations,  it  is  extremely  important  to 
make  sure  that  the  structural  properties  are  accurately  modeled  using  the  force  field.  As  shown 
in  Tables  1,  2,  and  3,  the  present  force  field  yields  excellent  agreement  with  experiment  and 
high-level  ab  initio  calculations  on  methyl  and  ethyl  nitrate.  The  average  deviation  of  the 
experimental  bond  lengths,  obtained  using  the  force  field,  is  0.5%  for  methyl  nitrate  and  0.6% 
for  each  of  the  conformers  of  ethyl  nitrate.  The  bond  angles  differ  on  average  by  1.4%  for 
methyl  nitrate,  1.5%  for  the  gauche  conformer  of  ethyl  nitrate,  and  1.3%  for  the  trans  conformer. 
For  the  larger  molecules  (NG  and  DEGDN),  similar  results  were  obtained.  A  graph  of  the  force- 
field-calculated  bond  lengths  and  bond  angles  of  DEGDN  and  NG  are  plotted  vs.  the  B3LYP/6- 
31  lg(d,p)  results  in  Figures  2  and  3,  respectively.  There  are  a  total  of  39  data  points  for  the  bond 
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Figure  2.  Comparison  of  the  Bond  Lengths  (in  A)  of  NG  and  DEGDN  Calculated  Using 
the  Force  Field  and  the  B3LYP/6-311g(d,p)  Reference. 


lengths  and  21  for  the  bond  angles.  As  indicated  in  the  figures,  excellent  agreement  between  the 
calculated  (force  field)  and  reference  (B3LYP)  data  is  obtained.  The  average  percentage 
deviation  in  the  bond  lengths  is  0.8%,  and  the  maximum  percentage  deviation  is  2.1%.  The  root 
mean  squares  (RMS)  percentage  deviation  is  0.9%.  For  the  bond  angles,  the  average  percentage 
deviation  is  1.1%,  with  a  maximum  percentage  deviation  of  3.0%.  The  RMS  percentage 
deviation  is  1.3%.  These  results  are  consistent  with  the  results  obtained  on  other  molecules 
using  the  COMPASS  force  field  [1, 2]. 

4.1.2  Vibrational  Frequencies.  A  molecular  mechanics  (MM)  force  field  is  different  firom  a 
spectroscopic  force  field.  Generally  speaking,  by  simultaneously  fitting  various  properties 
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Figure  3.  Comparison  of  the  Bond  Angles  (in  Degrees)  of  NG  and  DEGDN  Calculated 
Using  the  Force  Field  and  the  B3LinP/6-311g(d,p)  Reference. 

including  structures,  conformations,  and  vibrational  frequencies,  a  MM  force  field  can  only 
predict  the  vibrational  frequencies  with  a  modest  level  of  accuracy  (a  RMS  deviation  of 
approximately  20-50  cm'^).  Frequency  predictions  obtained  by  COMPASS  generally  fall  into 
this  category.  In  Tables  7-9,  the  vibrational  frequencies  of  methyl  and  ethyl  nitrate,  calculated 
using  the  force  field,  are  compared  with  experimentally  determined  frequencies  [20,  38]  and  with 
frequencies  computed  at  the  B3LYP/6-311g(d,p)  and  MP2/6-311g(d,p)  levels  of  theory.  In 
addition,  approximate  descriptions  of  the  modes  are  also  presented.  As  shown,  the  COMPASS- 
calculated  frequencies  agree  reasonably  well  with  the  experimental  values.  The  methyl  nitrate 
frequencies  have  an  average  deviation  of  -3.4  cm*'  from  the  experimental  values,  with  a 
maximum  deviation  of  94  cm'*.  The  RMS  deviation  in  the  methyl  nitrate  frequencies  is 
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Table  7.  Calculated  and  Experimental  Vibrational  Frequencies  (cm'^)  of  Methyl  Nitrate 


Mode 

Experimental* 

MP2/6-311g(d,p) 

B3LYP/6-311g(d,p) 

COMPASS 

Approx. 

Description* 

1 

3034 

3223 

3153 

2988 

CHj  stretch 

2 

3015 

3196 

3130 

2987 

CH3  stretch 

3 

2966 

3099 

3049 

2902 

CH3  stretch 

4 

1668 

1877 

1741 

1727 

NO2  stretch 

5 

1469 

1532 

1501 

1469 

CH3  deformation 

6 

1457 

1497 

1472 

1453 

CH3  deformation 

7 

1435 

1491 

1463 

1435 

CH3  deformation 

8 

1291 

1323 

1333 

1385 

NO2  stretch 

9 

1178 

1213 

1190 

1151 

CH3rock 

10 

1156 

1198 

1167 

1133 

CH3  rock 

11 

1019 

1081 

1020 

1113 

CO  stretch 

12 

855 

875 

871 

944 

ON  stretch 

13 

761 

767 

769 

743 

NO2  rock 

14 

658 

688 

663 

604 

NO2  deformation 

15 

571 

592 

572 

484 

NO2  wag 

16 

344 

357 

340 

310 

CON  bend 

17 

204 

234 

200 

212 

CH3  torsion 

18 

134 

138 

141 

113 

NO2  torsion 

‘Reference  [38]. 


52.6  cm'\  For  the  ethyl  nitrate  (including  both  the  trans  and  gauche  conformers)  frequencies,  the 
average  deviation  is  —13.2  cm'\  and  the  maximum  deviation  is  81  cm  The  RMS  deviation  is 
40.1  cm*'.  The  statistical  analysis  was  computed  for  18  methyl  nitrate  frequencies  and  54  ethyl 
nitrate  frequencies.  For  NG  and  DEGDN,  the  force  field  and  B3LYP/6-311g(d,p)  results  were 
compared  due  to  the  lack  of  experimental  data.  In  Figure  4,  the  COMPASS  frequencies  are 
plotted  against  the  B3LYP  frequencies.  There  are  111  data  points  plotted  with  an  average 
deviation  of  —12.1  cm'\  The  maximum  absolute  deviation  is  -190.7  cm'\  and  the  RMS 
deviation  is  65.4  cm'\  The  largest  deviations  are  found  in  the  hi^  frequency  region 
(>3000  cm’'),  which  correspond  to  the  C-H  stretch  modes.  By  comparing  the 
B3LYP/6-31  lg(d,p)  frequencies  against  the  experimental  values  for  methyl  and  ethyl  nitrates  (as 
shown  in  Tables  7-9),  it  is  clear  that  these  deviations  are  largely  due  to  the  over  estimates  of  C-H 
frequencies  by  the  B3LYP/6-31  lg(d,p)  method. 
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Table  8.  Calculated  and  Experimental  Vibrational  Frequencies  (cm'*)  of  Ethyl  Nitrate 
(Gauche  Conformer) 


Mode 

Experimental® 

MP2/6-311g(d,p) 

B3LYP/6-311g(d,p) 

COMPASS 

Approx. 

Description® 

1 

2990 

3203 

3141 

2971 

CH2  stretch 

2 

2986 

3183 

3134 

2970 

CH3  stretch 

3 

2977 

3178 

3115 

2964 

CH2  stretch 

4 

2949 

3124 

3108 

2913 

CH3  stretch 

5 

2940 

3074 

2901 

CH3  stretch 

6 

1660 

1871 

1734 

1741 

NO2  stretch 

7 

1519 

1529 

1508 

1522 

CH2  deformation 

8 

1460 

1508 

1488 

1472 

CH3  deformation 

9 

1444 

1505 

1482 

1461 

CH3  deformation 

10 

1392 

1444 

1419 

1413 

CH3  deformation 

11 

1369 

1414 

1401 

1397 

CH2  wag 

12 

1298 

1347 

1328 

1368 

CH2  twist 

13 

1288 

1317 

1323 

1283 

NO2  stretch 

14 

1270 

1205 

1183 

1192 

CH2rock 

15 

1161 

1133 

1109 

1140 

CH3  rock 

16 

1092 

1091 

1029 

1045 

CC  stretch 

17 

1025 

939 

910 

969 

CO  stretch 

18 

903 

868 

860 

873 

CH2  rock 

19 

809 

819 

807 

775 

NO  stretch 

20 

765 

765 

767 

744 

21 

643 

765 

649 

612 

NO2  scissor 

22 

573 

677 

574 

500 

NO2  rock 

23 

410 

419 

409 

375 

CCO  bend 

24 

346 

364 

349 

330 

CON  bend 

25 

199 

229 

225 

221 

CH3  torsion 

26 

no 

142 

124 

132 

CO  torsion 

27 

96 

94 

97 

90 

NO2  torsion 

“Reference  [20]. 


4.1.3  Conformations.  The  equilibrium  confirmations  of  methyl  and  ethyl  nitrate  are 
illustrated  in  Figure  1 .  The  NO2  group  is  coplanar  with  the  C-O-N  plane.  The  N-0  bond  has  a 
partial  double  bond  character,  as  indicated  by  the  relatively  high  energy  barrier  of  rotation  for  the 
NO2  group  arormd  the  N-O  bond,  as  given  in  Figures  5  and  6.  For  both  methyl  and  ethyl  nitrate, 
the  energy  barrier  heists  are  about  7  kcal/mol.  This  barrier  is  much  lower  than  what  is  observed 
for  a  true  double  bond,  but  is  significantly  higher  than  that  of  a  single  bond.  As  seen  in  Figures  5 
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Table  9.  Calculated  and  Experimental  Vibrational  Frequencies  (cm'*)  of  Ethyl  Nitrate 
(Trans  Conformer) 


Mode 

Experimental* 

MP2/6-311g(d,p) 

B3LYP/6-311g(d,p) 

COMPASS 

Approx. 

Description* 

1 

2990 

3197 

3126 

2972 

CH2  stretch 

2 

2986 

3185 

3112 

2969 

CHs  stretch 

3 

2977 

3162 

3098 

2966 

CH2  stretch 

4 

2949 

3104 

3055 

2912 

CH3  stretch 

5 

2940 

3093 

3044 

2900 

CH3  stretch 

6 

1660 

1873 

1732 

1741 

NO2  stretch 

7 

1519 

1549 

1525 

1496 

CH2  deformation 

8 

1460 

1522 

1501 

1468 

CH3  deformation 

9 

1444 

1503 

1484 

1457 

CH3  deformation 

10 

1392 

1448 

1426 

1410 

CH3  deformation 

11 

1369 

1413 

1403 

1397 

CH2  wag 

12 

1298 

1319 

1324 

1369 

CH2  twist 

13 

1288 

1316 

1294 

1299 

NO2  stretch 

14 

1270 

1205 

1180 

1197 

CH2  rock 

15 

1161 

1170 

1142 

1125 

CH3  rock 

16 

1122 

1085  ^ 

1028 

1063 

CC  stretch 

17 

1025 

954 

923 

990 

CO  stretch 

18 

903 

877 

876 

883 

CH2  rock 

19 

851 

843 

829 

777 

NO  stretch 

20 

765 

766 

769 

738 

NO2  wag 

21 

702 

732 

707 

639 

NO2  scissor 

22 

566 

583 

569 

494 

NO2  rock 

23 

378 

386 

373 

333 

CCO  bend 

24 

242 

263 

254 

255 

CON  bend 

25 

203 

233 

225 

215 

CH3  torsion 

26 

120 

125 

130 

118 

NO2  torsion 

27 

112 

101 

88 

92 

CO  torsion 

“Reference  [20]. 


and  6,  the  COMPASS  energy  profile  is  in  excellent  agreement  with  the  MP2/6-311g(d,p)  and 
B3LYP/6-31  lg(d,p)  results. 

Figures  7  and  8  are  the  energy  profiles  for  file  rotations  about  the  C-0  bond  in  methyl  and 
ethyl  nitrate,  respectively.  The  profile  of  methyl  nitrate  shows  a  low  energy  bamer  height  of 
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Reference  Frequencies  [B3LYP/6-311g(d,p)] 

Figure  4.  Comparison  of  the  Vibrational  Frequencies  (in  cm‘^)  of  NG  and  DEGDN 
Calculated  Using  the  Force  Field  and  the  B3LYP/6-311  g(d,p)  Reference. 

about  2.5  kcal/mol,  which  is  a  typical  value  for  a  single  bond.  The  rotation  about  the  same  C-0 
bond  in  ethyl  nitrate,  however,  has  a  more  complicated  pattern.  The  global  minimum  structure 
has  a  C-C-O-N  dihedral  angle  of  1 80°,  which  corresponds  to  the  trans  geometry.  In  addition, 
there  is  a  local  minimum  located  at  a  C-C-O-N  dihedral  angle  of  about  80°  (gauche  position). 
The  trans  to  gauche  isomerization  energy  barrier  is  approximately  2.0  kcal/mol;  the  transition 
state  to  this  isomerization  reaction  has  a  C-C-O-N  dihedral  angle  of  roughly  120°.  The  guache  to 
gauche  isomerization  has  a  barrier  height  of  about  10  kcal/mol,  and  the  transition  state  to  this 
isomerization  has  a  C-C-O-N  dihedral  angle  of  0°.  These  results  are  in  excellent  agreement  with 
the  HF/6-31g(d)  calculations  of  Durig  and  Sheehan  [20].  Apparently,  the  high  energy  is  largely 
due  to  nonbonded  interactions  (repulsion)  between  the  methyl  and  NO2  groups. 
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Figure  5.  Energy  (kcal/mol)  Profile  of  the  Rotation  of  the  NO2  Group  About  the  O-N 
Bond  of  Methyl  Nitrate. 

4.2  Liquid  Properties.  With  the  final  parameters,  NPT  simulations  were  performed  to 
calculate  the  densities  of  the  five  selected  liquid  nitrate  compounds  at  room  temperature  and  at 
temperatures  at  or  near  their  boiling  point.  The  boiling  point  temperature  for  NG  is  not  available 
in  the  literature;  however,  Stull  [39]  has  measured  the  vapor  pressure  of  NG  over  a  range  of 
temperatures  fi’om  400—524  K.  This  implies  that  NG  is  stable  at  these  temperatures  during  the 
time  fiume  in  which  the  measurements  were  made.  We  arbitrarily  chose  a  temperature  of 450  K 
to  run  the  elevated-temperature  NG  simulation,  which  is  approximately  in  the  middle  of  the 
range  over  which  the  experiment  was  conducted;  the  results  are  listed  in  Table  10.  The 
experimental  densities  measured  at  room  temperature  are  available  for  direct  comparison.  As 


22 


10.0 


Figure  6.  Energy  (kcal/mol)  Profile  of  the  Rotation  of  the  NOj  Group  About  the  0-N  Bond 
of  Ethyl  Nitrate. 

shown  in  the  table,  the  agreements  are  excellent  for  most  systems.  The  3%  deviation  for 
isopropyl  nitrate  is  slightly  larger  than  normal.  The  reason  is  unclear  at  this  point.  We  have  not 
foxmd  any  density  information  at  the  boiling  point  for  direct  comparison. 

The  values  of  the  heats  of  vaporization  were  calculated  using  equation  5,  with  simulated 
densities  at  two  temperatures  for  each  of  the  molecular  liquids.  Most  of  the  experimental  data 
are  derived  from  vapor  pressures  measured  at  room  temperature.  Our  results  calculated  at  room 
temperature  and  at  elevated  temperatures  are  in  excellent  agreement  with  the  experimental  data. 
The  deviations  are  within  the  same  normal  range  foxmd  in  many  molecules  that  COMPASS 
covers. 
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Figure  7.  Energy  (kcal/mol)  Profile  of  the  Rotation  of  the  CH3  Group  About  the  C-O  Bond 
of  Methyl  Nitrate. 

4.3  Crystal  Structure  of  Pentaerythritol  Tetranitrate  (PETN).  Although  the  nonbonded 
parameters  were  optimized  using  liquid  data,  it  has  been  shown  that  the  COMPASS  force  field  is 
capable  of  predicting  the  thermophysical  properties  of  both  liquids  and  crystals  over  an  extensive 
range  of  experimental  conditions.  For  additional  validation,  the  present  force  field  was  used  to 
calculate  the  crystal  structure  of  PETN.  The  unit  cell  was  constructed  from  the  experimental 
data  of  Trotter  [40].  The  unit  cell  contains  two  independent  molecules,  with  space-group 

J^=l^c  Figure  9  illustrates  the  projection  of  the  unit  cell  on  the  a-b  plane.  For  comparison 
purposes,  we  carried  out  both  energy  minimizations  and  NPT  simulations.  To  rigorously  test  the 
force  field  parameters,  the  calculations  were  performed  without  any  symmetry  constraints  (PI 
space-group).  In  addition,  all  of  the  cell  dimensions  and  angle  parameters  were  fully  relaxed. 
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Figure  8.  Energy  (kcal/mol)  Profile  of  the  Rotation  of  the  CH3CH2  Group  About  the  C-O 
Bond  of  Ethyl  Nitrate. 

The  MD  simulations  were  done  on  a  super  cell  that  contained  2x2x3  umt  cells.  The  energy 
minimization  was  performed  on  the  unit  cell.  In  all  simulations,  the  Ewald  summation  method 
was  used  for  both  the  electrostatic  and  the  vdW  energies. 

The  simulated  cell  parameters,  densities,  and  lattice  energies  are  listed  in  Table  11,  together 
with  the  experimental  data  of  Trotter  [40].  As  reported  before,  the  energy  minimization  method, 
which  corresponds  to  a  temperature  of  0°  K,  overestimates  the  density  by  about  5%.  The 
computed  density  from  a  room  temperature  NPT  simulation  agrees  well  with  the  experimental 
data.  The  calculated  sublimation  energies  are  also  given  in  this  table.  The  experimental  data  are 
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At  Room  Temperature 
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DEGDN  196.12  430.2 


Figure  9.  Projection  of  the  Unit  Cell  of  the  PETN  Crystal  on  the  a-b  Plane.  In  This  Figure, 
Red  Corresponds  to  Oxygen  Atoms,  Blue  Corresponds  to  Nitrogen  Atoms,  Black 
Corresponds  to  Carbon  Atoms,  and  White  Corresponds  to  Hydrogen  Atoms. 

reported  for  the  temperature  range  of  328-405  K;  our  NPT  simulation  results  (293  K)  agree 
reasonably  well  with  the  experimental  value. 

In  addition,  elastic  constants  of  the  PETN  crystal  were  calculated  using  the  analytical  second 
derivative  method.  With  the  energy-minimized  unit  cell,  the  calculated  values  are  significantly 
larger  than  the  experimental  data,  as  shown  in  Table  11.  Considering  that  the 
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Table  11.  Cell  Parameters,  Densities,  Sublimation  Energies,  and  Elastic  Constants  of 
PETN 


Energy  Minimization  (0  K) 

NPT  Dynamics  (293  K) 

Experiment® 

Unit  Cell  Length  (A) 

a 

9.117 

9.35 

9.38 

b 

9.117 

9.35 

9.38 

c 

6.721 

6.67 

6.69 

Unit  Cell  Angle  O 

a 

90 

90 

90 

90 

90 

90 

7 

90 

90 

90 

Density  (g/cm^) 

1.8811 

1.7997 

1.7837 

Sublimation  Energy  (kcal/mol) 

43.58 

40.49 

35.95 

(328-405  K) 

Elastic  Constants  (GPa) 

Energy  Minimized  Cell 

Cell  Dimensions  Fixed 

Experiment® 

cll 

21.6 

13.2 

17 

cl2 

11.3 

7.6 

5.4 

c33 

16.6 

11.1 

n/a 

c44 

8.9 

6.4 

5 

c66 

6.6 

6.0 

n/a 

“Experiment  from  reference  [40]. 


calculation  is  done  on  a  cell  in  which  the  density  has  been  overestimated,  this  discrepancy  is 
understandable.  To  mimic  the  experimental  conditions,  we  calculated  the  elastic  matrix  with  cell 
dimensions  and  angles  fixed  at  the  experimental  value.  Using  this  eonstraint,  the  results  agree 
well  with  the  measurement. 


5.  Conclusions 


The  COMPASS  force  field  has  been  extended  to  inelude  nitrate  esters  by  combining  the 
results  of  ab  initio  calculations  and  empirical  fitting  of  the  condensed-phase  properties  of  methyl 
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and  ethyl  nitrate.  The  ab  initio  derived  partial  charges  and  geometrical  parameters  (bond 
lengths,  bond  angles,  dihedral  angles,  and  torsion  angles)  were  parameterized  to  generate  the 
valence  parameters  in  the  force  field.  The  force  field  van  der  Waals  terms  were  optimized  using 
MD  simulations  of  the  thermophysical  properties  of  liquid  methyl  and  ethyl  nitrate.  Using  the 
resulting  force  field,  molecular  mechanics  calculations  have  been  performed  on  isolated  nitrates. 
The  calculated  molecular  structures  and  vibrational  frequencies  agree  well  with  experiment, 
high-level  ab  initio,  and  DFT  results.  In  addition,  room  temperatme  simulations  of  the  liquid 
properties  of  methyl  nitrate,  ethyl  nitrate,  propyl  nitrate,  butyl  nitrate,  isopropyl  nitrate, 
nitroglycerin,  and  diethyleneglycol  dinitrate  have  been  completed.  The  results  of  these 
simulations  agree  well  with  the  experimental  data.  Finally,  the  crystal  structure,  density, 
sublimation  energy,  and  elastic  constants  of  PETN  have  been  calculated  and  compared  to  the 
experiment.  These  results  suggest  that  this  new  force  field  is  capable  of  predicting  a  range  of 
properties  of  molecules  containing  the  nitrate  ester  functionality  in  both  condensed  and  gas 
phases. 
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New  Atom  Types  and  Parameters  for  Nitro  Compounds  and  Nitrates 

Atom  types 


Type 

Mass 

Element 

Connection 

n3o 

14.00674 

N 

3 

0I2 

15.99940 

0 

1 

o2n 

15.99940 

0 

2 

Equivalences 

Tvoe  NonB 

Bond 

An&le 

Torsion 

OOP 

n3o 

n3o 

n3o 

n3o 

n3o 

n3o 

0I2 

0I2 

oi  = 

ol= 

ol= 

ol= 

o2n 

o2n 

o2n 

o2n 

o2 

o2 

Bond  Increments 

c3a 

n3o 

0.2390 

-0.2390 

c4 

n3o 

0.2100 

-0.2100 

c4 

o2n 

0.3170 

-0.3170 

c4o 

n3o 

0.2100 

-0.2100 

c4o 

o2n 

0.3170 

-0.3170 

hi 

n3o 

0.1880 

-0.1880 

n3o 

ol= 

0.4280 

-0.4280 

n3o 

o2n 

0.0010 

-0.0010 

Bond 

I 

J 

RO 

M 

K3 

K4 

c3a 

n3o 

1.4300 

313.8329 

-568.6087 

600.9597 

c4 

n3o 

1.4740 

301.6051 

-535.7028 

555.0420 

c4 

o2n 

1.4350 

400.3954 

-835.1951 

1313.0142 

c4o 

n3o 

1.4740 

301.6051 

-535.7028 

555.0420 

c4o 

o2n 

1.4350 

400.3954 

-835.1951 

1313.0142 

hi 

n3o 

1.0400 

439.9346 

-943.7307 

1180.9318 

n3o 

ol= 

1.2100 

765.0664 

-2070.2830 

2793.3218 

n3o 

o2n 

1.4020 

300.0000 

-1000.0000 

2000.0000 

Comment 

nitrogen  in  nitro  group 
oxygen  in  nitro  group  (-NO2) 
oxygen  in  nitrates 
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Angle 


I 

J 

K 

ThetaO 

^  IG 

c3a 

c3a 

n3o 

118.8000 

29.2436  -8.8495  -6.6020 

hi 

c4 

n3o 

107.0000 

54.9318  -9.1333  -11.5434 

hi 

c4 

o2n 

108.7280 

58.5446  -10.8088  -12.4006 

c3a 

n3o 

ol= 

117.7000 

63.9404  -18.4524  -14.3129 

c4 

n3o 

ol= 

117.5000 

64.5228  -18.4582  -14.4215 

hi 

n3o 

ol= 

115.7000 

53.8034  -14.1991  -11.8708 

ol= 

n3o 

ol= 

128.0000 

95.1035  -47.4240  -27.9164 

c4 

o2n 

n3o 

108.5000 

55.7454  -10.0067  -6.2729 

c4 

c4 

o2n 

105.0000 

54.5381  -8.3642  -13.0838 

o2n 

n3o 

ol= 

112.8000 

85.5228  -18.4582  -14.4215 

Torsion 

I 

J 

K 

L 

vri) 

V(2) 

Ym 

c3a 

c3a 

c3a 

n3o 

0.0000 

7.2124 

0.0000 

hi 

c3a 

c3a 

n3o 

0.0000 

2.9126 

0.0000 

c3a 

c3a 

n3o 

ol= 

0.0000 

1.1600 

0.0000 

c4 

c4 

n3o 

ol- 

0.0000 

0.0000 

-0.3500 

hi 

c4 

n3o 

ol= 

0.0000 

0.0000 

-0.3500 

c4 

c4 

o2 

n3o 

0.0000 

-0.4000 

-0.2000 

ol= 

n3o 

o2 

c4 

0.0000 

2.0000 

0.0000 

Wilson  out  of  Diane 

I 

J 

K 

L 

KChi 

ChiO 

c3a 

c3a 

c3a 

n3o 

0.9194 

0.0000 

c3a 

n3o 

ol= 

ol= 

36.2612 

0.0000 

c4 

n3o 

ol= 

ol= 

44.3062 

0.0000 

hi 

n3o 

ol= 

ol= 

38.5581 

0.0000 

ol= 

n3o 

ol= 

o2 

45.0000 

0.0000 

Nonbond(X,J9-6> 


I 

r 

8 

n3o 

3.7600 

0.04800 

ol2 

3.4000 

0.04800 

o2n 

3.6500 

0.20000 
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Bond-bond 


I 

J 

K 

Ktb.b’) 

c3a 

c3a 

n3o 

21.0495 

c4 

c4 

o2n 

11.4318 

hi 

c4 

n3o 

3.3770 

hi 

c4 

o2n 

23.1979 

c3a 

n3o 

ol= 

93.7948 

o2n 

n3o 

ol= 

80.0000 

c4 

n3o 

ol= 

48.1403 

hi 

n3o 

ol= 

14.8266 

ol= 

ii3o 

ol= 

265.7106 

Bond-angle 


I 

J 

K 

Kfb^heta) 

Ktb'.theta) 

c3a 

c3a 

ii3o 

30.5211 

59.8025 

c4 

c4 

o2n 

2.6868 

20.4033 

hi 

c4 

ii3o 

12.2491 

30.5314 

hi 

c4 

o2n 

4.6189 

55.3270 

c3a 

n3o 

ol= 

40.3757 

92.1955 

c4 
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